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FOREWORD 


This  report  has  been  prepared  to  present  a  formal  description  of  the  Sun-Moon- 
Transformation  angle  support  software  package.  This  package  constructs  a  magnetic 
tape  file  of  inertial  solar  and  lunar  position  vectors  and  nutation  angles.  Values 
for  the  obliquity  of  the  ecliptic  and  the  equation  of  the  equinoxes  are  also 
computed  and  stored.  These  data  may  be  generated  in  either  the  B1 950.0  or  the 
J2000.0  reference  systems  and  are  intended  for  use  by  the  PULSAR  Doppler  satellite 
tracking  data  editor. 
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SECTION  1 
INTRODUCTION 

This  document  provides  a  basic  overview  of  the  computational  functions  of 
the  sun-moon-transformation  angle  support  software  (SMTS).  This  software 
package  was  designed,  documented  and  implemented  by  members  of  the  Space 
Flight  Sciences  Branch  of  the  Space  and  Surface  Systems  Division  at  the 
Naval  Surface  Weapons  Center  (NSWC)/Dahlgren,  Virginia. 

The  primary  function  of  the  SMTS  package  is  to  construct  a  tape  file  con¬ 
taining  inertial  solar  (r„,,  )  and  lunar  (r  )  position  vectors  and  the 

sun  moon 

transformation  angles  ( Aip ,  e,  Ac,  AH)  needed  to  rotate  between  the  earth- 
centered  inertial  and  the  earth-centered  earth-fixed  coordinate  reference 
frames.  These  data  may  be  generated  in  either  the  B1 950.0  or  the  J2000.0 
reference  systems.  The  transformation  angles  and  the  solar  position  vectors 

i- 

are  generated  once  a  day  at  o  ,  whereas  the  lunar  position  vectors  are  gener- 
ated  twice  a  day  at  on  and  12  .  The  resulting  tape  file  is  intended  for  use 
by  the  PULSAR  Doppler  Satellite  tracking  data  editor1. 

A  functional  overview  of  the  SMTS  package  is  shown  in  Figure  1-1.  As  can  be 
seen  from  this  figure,  the  SMTS  package  is  comprised  of  three  basic  computa¬ 
tional  functions:  the  process  control  function  (PCF),  the  sun-moon  ephemeris 
generator  (SMEG),  and  the  transformation  angle  generator  (TAG).  Also  shown 
on  this  figure  is  the  primary  inter-function  data  flow.  A  descriptive  over¬ 
view  of  each  of  these  functions  is  presented  in  the  following  sections. 

As  input  the  SMTS  package  requires  a  standard  Jet  Propulsion  Laboratory  (JPL) 
export  tape.  Also  needed  are  the  year  (YRB)  and  day  number  (DAYB)  to  begin 
the  sun-moon  file;  the  year  (YRE)  and  day  number  (DAYE)  to  end  the  sun-moon 


1  Parks,  A.D.  and  Hicks,  T.I.,  A  Mathematical  Description  of  the  PULSAR 
Doppler  Satellite  Tracking  Data  Editor,  NSWC/DL  TR  82-391  , 
Dahlgren,  Virginia,  1982. 
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1.0  INTRODUCTION  (con't) 

file;  a  switch  (IREF)  to  select  either  the  J2000.0  or  B1950.0  reference 
system;  and  a  flag  ( I JPL )  to  inform  SMTS  which  reference  system  (i.e.  J2000.0 
or  B1 950.0)  the  JPL  export  tape  is  based  upon. 

In  order  to  initiate  its  processing,  the  SMTS  reads  a  single  data  card 
specifying  YRB,  DAYS,  YRE,  DAYE,  IREF,  and  IJPL.  The  JPL  tape  is  accessed 
and  sun  and  moon  position  vectors  are  generated  for  the  entire  time  span 
using  Chebyshev  expansions.  These  vectors  are  then  stored  in  a  temporary 
file  (TAPE  1 )  .  The  transformation  angles  are  generated  on  a  daily  basis 
and  merged  with  the  associated  sun  and  moon  position  vectors  on  the  final 
sun-moon-transformation  angle  tape  (SMT). 
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SECTION  2 

THE  PROCESS  CONTROL  FUNCTION  (PCF) 

2.1  FUNCTIONAL  DESCRIPTION 

The  principal  tasks  performed  by  the  PCF  are  to  receive  and  retrieve  data, 
direct  processing  flow,  and  output  the  desired  SMT  file.  Specifically  the 
PCF: 

U)  receives  data  from  the  input  card; 

[Li)  converts  input  times  to  Julian  days; 

{Lit)  determines  the  time  span  for  which  to  construct  a  SMT  tape; 

(.tv)  calls  SMEG  to  produce  sun-moon  position  vectors  for  the  re¬ 
quired  time  span  and  write  them  to  a  temporary  file; 

(v)  calls  TAG  to  compute  transformation  angles  for  the  span  of 
interest;  and 

(u^)  merge  the  sun-moon  position  vectors  with  the  transformation 
angles  to  form  the  SMT  tape. 

The  flow  of  the  PCF  function  is  presented  in  Figure  2-1. 

2.2  PROCESSING  EQUATIONS  AND  SMT  FILE  FORMAT 

The  input  data  is  read  from  the  single  input  card  and  the  beginning  Julian 
day,  JD3,  and  processing  span,  AT,  are  computed  as  follows: 

JDg  =  367.  YRB  -  INT  [(7  YRB)/4.]  +  DAYB  +  1721043.5  (2.1) 

and 


AT  =  366 .(YRE  -  YRB  +  1 )  .  (2.2) 

The  time  in  Julian  Centuries  used  by  the  TAG  function  to  generate  the  asso¬ 
ciated  transformation  angles  is  computed  using 
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FIGURE  2-1.  PROCESS  CONTROL  FUNCTION  LOGIC  FLOW 


TR  82-389 


2.2  PROCESSING  EQUATIONS  AND  SMT  FILE  FORMAT  (con ' t) 


T 


JDc  -  2451545.0 


36525. 


(2.3) 


where  JDC  is  the  current  Julian  Ephemeris  Date  for  the  time  of  interest  and 
is  computed  using  an  expression  analagous  to  equation  (2.1).  This  is  the 
same  t  used  in  SECTION  4.0. 


The  format  for  the  SMT  tape  is  given  in  Table  2-1.  It  should  be  mentioned 
that  the  pole  wander  corrections,  Ap  and  Aq,  and  the  UT1-UTC,  At,  correc¬ 
tions  are  zeroed  out  on  the  file.  This  is  done  because  these  quantities 
are  computed  and  predicted  internally  by  PULSAR  (see  SECTION  2.2.3  in 
reference  1 ) . 
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TABLE  4-2.  INTEGER  MULTIPLES  AND  COEFFICIENTS  REQUIRED  TO  COMPUTE 
Aty  AND  Ae  IN  THE  BT 950.0  REFERENCE  SYSTEM  (con't) 
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TABLE  4-2.  INTEGER  MULTIPLES  AND  COEFFICIENTS  REQUIRED  TO  COMPUTE 
Ai|/  AND  Ae  IN  THE  B1950.0  REFERENCE  SYSTEM  (con't) 


V 

X. 
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■t 

n3. 

X, 

n4. 

X- 

0  5  . 
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Ae 

3 

h 
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2 
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1 
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1 
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INTEGER  MULTIPLES  AND  COEFFICIENTS  REQUIRED  TO  COMPUTE 
A\l>  AND  Ae  IN  THE  B1950.0  REFERENCE  SYSTEM 
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h 
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TABLE  4-1.  CONSTANTS  REQUIRED  TO  COMPUTE  FUNDAMENTAL  SUN-MQQN  ANGLES  (*) 


C, 

c2 

.0518 

33.22454 

.064 

31.31 

-.0108 

-.5724 

-.012 
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0.0 

7.504 
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4.2  PROCESSING  EQUATIONS  (con't) 


e  =  e  +  Ae 


(4.8) 


where 


f  23.439280808+  [(5.036  X10_7t-  1.31  X10_7)t-  .013014269]t  ,  B1950.0 
e(°)  =|  (4.9) 


23.4392911111  +  [(5.036  X  10-7t  -  1 .639  X  10_7)t  -  .0130041667]t  ,  J2000.0 


and 


AH  =  tan"^  [tanAijj  cos  e] 


(4.10) 


-16- 


-A 


«  •  •  •  . 


TR  82-389 


4.2  PROCESSING  EQUATIONS  (con't) 


In  the  last  two  equations 


t  =  Julian  Ephemeris  Date,  J2000.0 


t  =  2451545.0 


,  J2000.0 


the  ny.  (/  =  1,  2,  ...,  5)  are  integers  and: 


l  is  the  mean  elongation  of  the  moon  from  the  sun; 


l  is  the  mean  argument  of  latitude  of  the  moon; 


F  is  the  mean  anomaly  of  the  sun; 


d  is  the  mean  anomaly  of  the  moon;  and 


ft  is  the  mean  longitude  of  the  ascending  node  of  the  moon. 

These  angular  quantities  are  computed  in  units  of  seconds  of  arc  from  expres 
sions  having  the  following  form: 


X  =  [ (Cjt  +  C2)t  +  C3]t  +  C4  +  1296000.000  mod  (CgT,  1.000). 


where  t  is  given  by  equation  (4.3)  and  the  constants  C i  {l  =  1,  2,  ...,  5) 
are  given  in  Table  4-1.  The  plus  sign  after  C4  applies  for  all  cases  except 
when  X  =  ft,  in  which  case  the  minus  sign  applies.  Values  for  the  a^> 
and  n;  of  equations  (4.1),  (4.2),  and  (4.4)  are  given  in  Tables  4-2  and  4-3 
for  the  B1950.0  and  J2000.0  reference  systems,  respectively. 


The  obliquity  of  the  ecliptic,  e,  and  the  equation  of  the  equinoxes,  AH, 
are  computed  using 
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FIGURE  4-1. 
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SECTION  4 

THE  TRANSFORMATION  ANGLE  GENERATOR  (TAG) 


4.1  FUNCTIONAL  DESCRIPTION 


The  transformation  angle  generator  computes  on  a  daily  basis  at  o*1  the 
nutation  in  longitude,  Aip;  the  nutation  in  obliquity,  Ae;  the  obliquity  of 
the  ecliptic,  e;  and  the  equation  of  the  equinoxes,  AH.  As  the  angles  are 


computed,  they  are  merged  with  rsun  and  ^moon  on  the  SMT  tape.  These  angles 


may  be  computed  in  either  the  B1950  or  the  J2000  reference  frame  by  setting 
the  IREF  flag.  The  process  flow  of  the  TAG  function  is  shown  in  Figure  4-1. 


4.2  PROCESSING  EQUATIONS 

The  nutation  in  longitude,  Aip,  and  the  nutation  in  obliquity,  Ae,  in 
seconds  of  arc  are  computed  in  both  reference  systems  using  the  following 
expansions: 


N 


Aip  =  10  (<x4  +  B^t)  sin 

1=  1 


0£ 


(4.1) 


and 


Ae  =  10-4£(y^)  cos  0^.  , 

.0=1 


(4.2) 


(N  =  69  for  the  B1950.0  system  and  N  =  106  for  the  J2000.0  system)  where 


t  -  t. 


T  = 


36525.0 


(4.3) 


and 


n  Z  +  n2  l  +  n3  F  +  n4  d  +  n5  ft 

A*  A,  A*  A_ 


(4.4) 
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3.2  PROCESSING  EQUATIONS  (con't) 

where  the  &£yfe  and  are  the  Chebyshev  expansion  coefficients.  The  velocity 
and  acceleration  polynomials  are  obtained  by  differentiating  recursion  relation 
(3.4),  i.e. 


Tn+1  "  2tc  "  ^n-1  +  2Tn 


(3.10) 


2t  T  -  T  !  +  4T 

c  n  n-1  n 


(3.11) 


As  was  mentioned  earlier,  if  desired,  it  is  possible  to  transform  the  sun  and 
moon  between  the  J2000.0  and  B1 950.0  reference  systems  by  setting  the  IJPL 
flag.  These  transformations  are  given  by2: 


r^(J2000)  =  R  r^(B1950)  ,  ( l  =  sun,  moon) 


(3.12) 


r^(B1950)  =  R  r^(J2000),  (•€  =  sun,  moon)  , 


(3.13) 


where 


0.9999256791774783  -  0.0111815116768724  -  0.0048590038154553 


R  =  0.0111815116939975  0.9999374845751042  -  0.0000271625775175  .  ( 


0.0048590037714450  -  0.0000271704492210  0.9999881946023742 


2  Standish,  Jr.,  E.  Myles,  Orientation  of  the  JPL  Eohemerides,  DE200/LE200, 


To  The  Dynamical  Equinox  of  J2000,  privately  obtained  pre-print,  April, 
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3.2  PROCESSING  EQUATIONS  (con't) 

Recursion  relations  are  used  to  generate  the  Chebyshev  polynomials  at  t 


Tn+1  =  2tc  Tn  "  Tn-1  * 


where 


1  , 


2  . 


The  j  position  components  for  the  sun  and  moon  at  t  are  given  by 


(rt)i 


X  am  T 


=  sun,  moon;  /  =  1 ,  2,  3) , 


where  is  the  associated  Chebyshev  expansion  coefficient  valid  for  the 
interval  [ts,  t£]  and  obtained  from  the  JPL  tape. 

It  should  be  mentioned  that,  although  they  are  not  used  to  construct  the 
SMT  tape,  the  velocity  and  acceleration  for  the  sun  and  moon  can  also  be 
obtained  from  the  SMEG  function: 

m 

(r£) .  =  y;  ife  ,  U  =  sun,  moon;  j  =  1,  2,  3)  1 


mi 

(rz)  ■  =  Yijk  Tfe  ,  U  =  sun,  moon;  j  =  1,  2,  3) 
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SECTION  3 

THE  SUN-MOON  EPHEMERIS  GENERATOR  (SMEG) 

3.1  FUNCTIONAL  DESIGN 

The  sun-moon  ephemeris  generation  function  is  us  d  by  the  SMTS  package  to 
provide  diurnal  solar  location  vectors  and  semi-diurnal  lunar  location  vec¬ 
tors.  These  vectors  are  computed  using  Chebyshev  expansions  and  are  stored 
on  a  temporary  data  file  before  being  merged  with  the  transformation  angle 
data.  It  should  be  mentioned  that  if  desired,  the  IJPL  flag  may  be  set  to 
transform  the  position  vectors  to  the  J2000.0  reference  system,  if  initially 
expressed  in  the  B1950.0  reference  system;  or  to  the  B1950.0  reference  system, 
if  initially  expressed  in  the  J2000.0  system.  Of  course,  care  must  be  taken 
to  insure  that  the  same  reference  system  is  being  used  consistently  for  the 
sun  and  moon  position  vectors  and  the  transformation  angles  when  constructing 
the  SMT  tape.  The  fundamental  SMEG  processing  flow  is  shown  in  Figure  3-1. 


3.2  PROCESSING  EQUATIONS 

Values  for  sun  and  moon  position  coordinates  are  obtained  in  SMEG  using  expan¬ 
sions  in  Chebyshev  polynomials.  The  coefficients  required  for  these  expansions 
are  read  from  the  JPL  tape  and  are  valid  over  certain  time  subintervals  of  the 
span  covered  by  the  tape.  For  a  time  t  in  a  subinterval  [ts,  t^],  i.e.  for 

t$  <  t  <  tE  ,  (3.1) 

a  Chebyshev  time,  tc,  is  computed  from 
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Table  2-1.  SMT  Tape  Format 


RECORD 

NUMBER  OF 
WORDS  IN 
RECORD 

VARIABLE 

DESCRIPTION 

Header 

3 

REC1 

Record  Type  =  110XXX*  for  the  B1950.0  system 

=  210XXX  for  the  J2000.0  system 

YRAt 

Year  Number 

yra2 

Number  of  days  in  block 

Data 

20 

REC2 

Record  Type  =  lllXXX  for  the  B1950.0  system 

=  211 XXX  for  the  J2000.0  system 

DAY 

Day  Number 

SUN 

L 

Coordinates  of  sun  at  0  (three  words) 

M00Nx 

Coordinates  of  moon  at  0h  (three  words) 

MOON  2 

Coordinates  of  moon  at  12  (three  words) 

DPSI 

Nutation  in  longitude  (radians) 

TRUEPS 

Obliquity  of  the  ecliptic  (radians) 

DEPS 

Nutation  in  obliquity  (radians) 

DTS 

UT1  -  UTC  timing  correction  (seconds) 

DELHT 

Equation  of  the  equinoxes  (radians) 

P 

Polar  displacement  toward  Greenwich  (radians) 

Q 

Polar  displacement  west  (radians) 

DUM  1 

Dummy  word 

DUM  2 

Dummy  word 

*  XXX  =  XXXth  block  of  taDP 


1 

V 

"u 

1 

0 

0 

2 

0 

0 

3 

0 

0 

4 

0 

0 

5 

0 

1 

6 

1 

0 

7 

0 

1 

8 

0 

0 

9 

1 

0 

10 

0 

-1 

11 

1 

0 

12 

0 

0 

13 

-1 

0 

14 

1 

0 

15 

0 

0 

16 

-1 

0 

17 

-1 

0 

18 

1 

0 

19 

2 

0 

20 

-2 

0 

21 

0 

0 

22 

2 

0 

23 

2 

0 

24 

1 

0 

0 

2 

2 

0 

0 

0 

2 

2 

2 

2 

0 

2 

2 

0 

0 

2 

0 

2 

0 

2 

2 

2 

0 

2 


0 

-2 

0 

0 

0 

0 

-2 

0 

0 

-2 

-2 

-2 

0 

0 

2 

2 

0 

0 

-2 

0 

2 

0 

0 

-2 


1 

2 

2 

2 

0 

0 

2 

1 

2 

2 

0 

1 

2 

1 

0 

2 

1 

1 

0 

1 

2 

2 


Ae 

*1 

U 

-171996. 

-174.2 

92025. 

8.9 

-13187. 

-1.6 

5736. 

-3.1 

-2274. 

-0.2 

977. 

-0.5 

2062. 

0.2 

-895. 

0.5 

1426. 

-3.4 

54. 

-0.1 

712. 

0.1 

-  7. 

0.0 

-517. 

1.2 

224. 

-0.6 

KO 

00 

CO 

i 

-0.4 

200. 

0.0 

-301. 

0.0 

129. 

-0.1 

217. 

-0.5 

-95. 

0.3 

-158. 

0.0 

-  1. 

0.0 

129. 

0.1 

-70. 

0.0 

123. 

0.0 

-53. 

0.0 

63. 

0.1 

-33. 

0.0 

63. 

0.0 

-  2. 

0.0 

-59. 

0.0 

26. 

0.0 

-58. 

-0.1 

32. 

0.0 

-51. 

0.0 

27. 

0.0 

48. 

0.0 

1. 

0.0 

46. 

0.0 

-24. 

0.0 

-38. 

0.0 

16. 

0.0 

-31. 

0.0 

13. 

0.0 

29. 

0.0 

-1. 

0.0 

29. 

0.0 

-12. 

0.0 

0 

0 
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TABLE  4-3.  INTEGER  MULTIPLES  AND  COEFFICIENTS  REQUIRED  TO  COMPUTE  A\p  AND 
Ae  IN  THE  J2000.0  REFERENCE  SYSTEM  (con’t) 
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TABLE  4-3.  INTEGER  MULTIPLES  AND  COEFFICIENTS  REQUIRED  TO  COMPUTE  Aip  AND 
Ae  IN  THE  J2000.0  REFERENCE  SYSTEM  (con't) 


n 

"u 

\ 

\ 

"u 

A\p 

Ae 

ai 

U 

73 

1 

0 

0 

2 

1 

-1. 

0.0 

0. 

0.0 

74 

1 

0 

2 

2 

1 

-1. 

0.0 

1. 

0.0 

75 

1 

1 

0 

-2 

1 

-1. 

0.0 

0. 

0.0 

76 

0 

1 

0 

2 

0 

-1. 

0.0 

0. 

0.0 

77 

0 

1 

2 

-2 

0 

-1. 

0.0 

0. 

0.0 

78 

0 

1 

-2 

2 

0 

-1. 

0.0 

0. 

0.0 

79 

1 

0 

-2 

2 

0 

-1. 

0.0 

0. 

0.0 

80 

1 

0 

-2 

-2 

0 

-1. 

0.0 

0. 

0.0 

81 

1 

0 

2 

-2 

0 

-1. 

0.0 

0. 

0.0 

82 

1 

0 

0 

g| 

0 

-1. 

0.0 

0. 

0.0 

83 

0 

0 

El 

0 

-1. 

0.0 

0. 

0.0 

84 

0 

2 

D 

2 

-1. 

0.0 

0. 

0.0 

85 

0 

2 

El 

2 

-1. 

0.0 

0. 

0.0 

86 

-2 

0 

2 

II 

2 

-1. 

0.0 

1. 

0.0 

87 

2 

0 

2 

H 

2 

-1. 

0.0 

0. 

0.0 

88 

0 

-1 

2 

0 

1 

-1. 

0.0 

0. 

0.0 

89 

0 

0 

-2 

0 

1 

-1. 

0.0 

0. 

0.0 

90 

0 

0 

IK 

-2 

2 

1. 

0.0 

0. 

0.0 

91 

0 

1 

I 

0 

2 

1. 

0.0 

0. 

0.0 

92 

1 

1 

2 

-2 

2 

1. 

0.0 

-1. 

0.0 

93 

0 

2 

-2 

2 

1. 

0.0 

0. 

0.0 

94 

0 

2 

2 

2 

1. 

0.0 

-1. 

0.0 

95 

-1 

0 

0 

0 

2 

1. 

0.0 

-1. 

0.0 

96 

0 

0 

-2 

2 

1 

1. 

0.0 

0. 

0.0 

TR  82-389 


TABLE  4-3.  INTEGER  MULTIPLES  AND  COEFFICIENTS  REQUIRED  TO  COMPUTE  Ai|>  AND 
Ae  IN  THE  J2000.0  REFERENCE  SYSTEM  (con't) 


r'T 


NSWC  TR  82-389 
DISTRIBUTION 


Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22314  (12) 

Library  of  Congress 

ATTN:  Gift  and  Exchange  Division  (4) 

Washington,  DC  20540 

National  Aeronautics  and  Space  Administration 

Scientific  and  Technical  Library 

Code  NHS  22,  Rm.  BA39 

600  Independence  Avenue,  SW 

Washington,  DC  20546  (2) 

Defense  Mapping  Agency 

ATTN:  Mr.  Jack  Calender  (10) 

Washington,  DC  20305 

Defense  Mapping  Agency 
Hydrographic/Topographic  Center 
ATTN:  Mr.  Patrick  Fell  (10) 

Washington,  DC  20390 

Defense  Mapping  Agency  Aerospace  Center 
ATTN:  Dr.  Robert  Ballew  (8) 

St.  Louis,  MO  63118 


Naval  Electronics  Systems  Command 
Navy  Space  Project,  PME106 

Washington,  DC  20360  (3) 

Office  of  Chief  of  Naval  Operations 
Naval  Oceanography  Division  (NOP-952) 

Bldg.  1,  U.S.  Naval  Observatory  (2) 

Washington,  DC  20390 

Office  of  Naval  Operations 

Navy  Space  Systems  Division  (NOP-943)  (2) 

Washington,  DC  20350 

Naval  Research  Laboratory 

ATTN:  Mr.  Al  Bartholomew  (3) 

Washington,  DC  20375 

Naval  Oceanographic  Office 

Bay  St.  Louis,  MS  39522  (2) 


(1  ) 


END 

FILMED 

9-85 


DTIC 


